lb = -5;
ub = 5;
n = 1000;
[psi,x] = gauswavf(lb,ub,n,8);
plot(x,psi)
title('Order 8 Gaussian Wavelet')
grid on